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ABSTRACT 



Aims. We present the results of the numerical simulations of the interaction between a magnetized star and an imperfectly conducting accretion 
disk. 

Methods. To analyze the "star-disk" interaction we numerically investigate the MHD equations used Godunov-type high resolution numerical 
method. 

Results. It was found that the "star - disk" interaction occurs with quasi-periodic reconnection of the magnetic field coronal loops and 
plasmoid ejections. In the case of the perfect disk conductivity the evolution of the coronal magnetic field leads to the periodic outflow of 
angular momentum from the disk. In the case of an imperfectly conducting disk the configuration of the magnetic field is formed such that the 
disk angular momentum carried by magnetic field gets balanced by angular momentum carried by matter. 

Key words. Magnetohydrodynamics (MHD) - accretion disks - stars: magnetic field - interstellar medium: evolution - methods: numerical 



1. Introduction 

This paper studies the evolution of the coronal magnetic field 
linked with a magnetized star and its accretion disk. We sug- 
gest that the plasma differential rotation along magnetic field 
lines is the reason of a "star-corona-disk" system evolution. 
We suppose that the magnetic field lines are frozen in the 
perfectly conducting coronal plasma. The differential rotation 
leads to the generation of the magnetic field toroidal com- 
ponent. Magnetic pressure increases in corona inner part and 
plasma is pushed towards the outer part together with the mag- 
netic field lines. As a result there is deformation or even open- 
ing of poloidal magnetic field lines adopting a new configura- 
tion. The type of this new configuration is determined by sev- 
eral factors. One of them is the electrical conductivity of the 
relatively cold disk plasma. 

We assume in this model that the imperfect plasma conduc- 
tivity is essential only in the disk and is determined by velocity 
turbulent fluctuations. The value of the turbulent magnetic dif- 
fusivity we consider as a free parameter of the problem. To ob- 
tain an acceptable range for this parameter, it is supposed that 
the coefficient of the turbulent magnetic diffusivity agrees in 
order of magnitude with the turbulent viscosity accepted in the 
standard ff-model of Shakura-Sunyaev accretion disk (Shakura 
& Sunvaev lT3751 . 



We consider that the disk is formed by relatively dense and 
cold matter. The disk is Keplerian and sound speed is much less 
than the Keplerian one. It means that disk is geometrically thin. 
In this model, the disk is considered as an infinitely thin, con- 
ductive plane. It should be mentioned that disk has a compli- 
cated structure and its interaction with the magnetic field does 
not reduce to magnetic compression and magnetic field lines 



slippage relatively to matter (Balbus et al. 1995). 

Besides the dynamics of the magnetic field, it is worthy to 
know its configuration after the opening of the magnetic field 
lines. This configuration defines the disk and the magnetic field 
evolution at large time scales. The main factor influencing this 
evolution is the rate of angular momentum transfer from the 
disk. The important role in this process belongs also to the mag- 
netic field (Feireira [T99"7l . 
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A sufficiently strong magnetic field can lead also to the for- 
mation of a matter outflow from the disk to the corona. In case 
of a thin Keplerian disk, a criterion for the "wind" formation 
beginning has been given by Blandford & Payne (1982i. To 
generate this outflow from a Keplerian disk, the magnetic field 
line should be inclined to the rotation axis with an angle of 
more than 30°. In that case a magnetic field line plays the role 
of a "rail" along which matter leaves the disk. Thus the final 
magnetic field configuration determines strongly both the rate 
of disk accretion, through the rate of angular momentum out- 
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flow, and the rate of matter outflow along the magnetic field 
lines inclined to the rotation axis. 

There are lots of papers concerning the evolution of the 
magnetic field in interaction with a disk. We would like to men- 
tion some of them. Hayashi et al. (1996 1 presented resistive 
MHD simulations with diffusive accretion disk and dipole stel- 
lar magnetic field topology. Kuwabara et al. (2000 1 introduced 
resistivity to simulate the effects of turbulent magnetic diffusiv- 



ity. Romanova et al. ( 2002 1 studied the disk accretion to a rotat- 
ing magnetized star with an aligned dipole moment and asso- 
ciated funnel flows, they described the "star-disk" interaction 
for the cases of fast and slowly rotating star. These simulations 
included a treatment of the disk vertical structure. Some pa- 
pers presented the results of ideal MHD simulations considered 
the disk like a boundary condition. Ustyugova et al. ( 1999jl for 
initially split-monopole magnetic field configuration, Ouyed & 
Pudritz ( 1997 1 for dipole topology, Krasnopolsky et al. (Q999) 
studied MHD outflows from the disk. Fendt & Elstner ( 1999 
2000 1 studied the "star-disk" interaction, they observed the 
opening of magnetic field lines and outflows. Fendt & Cemeljic 



( 2002 1 studied the jet formation and propagation by used the 
resistive MHD equations. 



Shu et al. (1994) presented steady state dynamics of ac- 
cretion from viscous and imperfectly conducting disk. Ferreira 
( 1997i) gave self-similar solution for the resistive disk for sta- 
tionary MHD equations. Uzdensky (2004) gave a review of 
modern theoretical scenario of the star-disk system interaction, 
including both stationary and time-dependent ones. 

In this paper we consider the evolution of the magnetic field 
configuration from an initial dipole-like topology into the final 
one. We analyze the time-dependent ideal MHD equations us- 
ing the high resolution Godunov-type method. We consider for 
a setup model of "star-disk" system that the generated flow 
is axisymmetric and symmetric in relation with the equatorial 
plane. The dependence of the disk angular momentum outflow 
rate on the disk surface electrical conductivity is investigated. 

The structure of the paper is as follows. After the intro- 
duction we propose the model and the evaluation of the disk 
surface magnetic diffusivity. Then we present the numerical 
method and describe our results for number disk surface mag- 
netic diffusivity. Finally we formulate the main conclusions. 



2. Statement of the problem 

Let us consider a "star-corona-disk" system governed by the 
magnetic field (see Fig.^. In our model the star has mass M», 
magnetic moment p* and angular velocity Q» . It is assumed that 
the star rotation axis aligns with its magnetic axis z. The disk 
rotating around the star is considered to be infinitely thin and 
located at z = 0. It is supposed that the disk is rotating follow- 
ing Keplerian orbits around the star and the disk is imperfectly 
conducting. The particle velocity in a Keplerian orbit with ra- 
dius r is Vk = y/GM* / r, where G - 0.667 ■ 10 _7 cm 3 g _I sec -2 is 
gravitational constant. Thus, the disk is differentially rotating. 



Q 
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Star R m R d 



Disk 



Rout R 



Fig. 1. Sketch of a magnetically linked "star-disk" system. The 
solid lines show the magnetic field lines. R m , R out are the edges 
of the computational domain. The star is inside the computa- 
tional domain, /?» < R m . 



The coronal plasma electro-conductivity is big enough, so 
we can describe the flow by the system of the ideal MHD equa- 
tions 



dp 

at 



+ V-(p«) = 0, 



^-+VT=pg, 

at 

— - V x (u x B) = 0, 

dt 



(1) 



dpS 



+ V-(pSH) = 0, 



dt 
VB = 0. 

Here T ik = p uiu k + p 5 ik + 



1 

4^ 



B 2 

-BiB k + —8 ik 



is the stress 



tensor; u is the plasma velocity; B is the magnetic field; p and 
p are the plasma density and pressure; S - pip 7 is the entropy 
function; y is the adiabatic index; g = -V(D g is the gravita- 
tional acceleration; cD g = -GM*/R is the star gravity potential; 
R is the distance from the gravitating center. 

The system ([0 is solved in spherical coordinates (R, <p, 8), 
with 8 being the polar angle with the symmetry axis. Velocity 
u and magnetic field B have all their components u = (u, v, w) 
and B = (B R , B v , Bg). 

The aim of this paper is to investigate the character of the 
magnetic field evolution and field topology depending on the 
disk surface magnetic diffusivity ( = c 2 link (A is surface elec- 
tric conductivity, c is the speed of light) at large time scales. As 
it will be shown below, it used in boundary conditions set on 
the equatorial plane z = 0. 
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2.1. Dimensionless variables and typical quantities for 
T Tauri stars 

The dimensionless form of Q is received in standard way. As 
distance scale, Rq, we take one third of the distance from the 
star center to the inner edge of the disk. Thus, in dimensionless 
units the disk inner edge is = 3R ln . The inner radius of the 
computational domain in dimensionless units is R[ n = 1 . The 
time and the velocity scales are chosen that in dimensionless 

units GM* = 1 . This requirement yields ?o = JR^/GM* as the 
time scale and vo = Ro/to as the velocity scale. As a magnetic 
field scale, Bo is taken and then density, pressure and magnetic 
moment scales are po = B^/v^, po = B 2 V p*o = BqRq. 

As typical quantities of T Tauri stars we adopt standard 
values like a star mass M» = 0.8Mq = 1.6 ■ 10 33 g, and 
Rd = 9Rq = 5.4 ■ 10 11 cm. Therefore the distance scale is 

= = 1.8 • 10 11 cm, the time and velocity scales are 

0.74 • 10 4 s and 2.43 • 10 7 cm s respectively. The Keplerian 
rotation period at the disk inner edge is 8.3 days. The simula- 
tion region size is 1.134 • 10 12 cm. 

The star magnetic moment is taken such that the magnetic 
field on the star surface is 300 G, so Bo = 6.5 G. Thus, on 
the disk inner part, the dipole magnetic field is 2.4 G, R t = 
3/5/?i n , p* = 10 and the magnetic moment is 3.8 ■ 10 35 Gcm 3 . 
The density scale is 1.44 • 10~ 14 gcrrr 3 , typical for the disks 
around T Tauri stars. 

2.2. Evaluation of the electric conductivity in the disk 

We suggest that turbulent diffusion of magnetic field is deter- 
mined by the same processes that determine turbulent viscosity, 
which leads to angular momentum transport in the disk. Thus, 
it is assumed that turbulent magnetic diffusivity T] t is in the or- 
der of turbulent viscosity like in the Shakura-Sunyaev model 
(Shakura & Sunyaev 1973): rj t = a t c s h, where c s is the sound 
speed in the disk, h is the disk half-thickness, a t is the dimen- 



sionless coefficient varying, according to Balbus (2003 1, in this 
range 0.01 4- 0.6. 

Under hydrostatic equilibrium, the Keplerian disk half- 
thickness h can be found from the relation: (h/r) 2 + b(h/r) - 
(c s /V k ) 2 = 0, where b = r(B 2 r + fi 2 )/(47r£V 2 ), S is the sur- 
face density, Vk is Keplerian velocity, B r , B^ are the magnetic 
field components (Bisnovatyi-Kogan & Lovelace 2001 1. In any 
case, even without taking into account the magnetic compres- 
sion, the disk half-thickness satisfies the inequality h S5 c s /D/;, 
where Qk is Keplerian angular velocity. Therefore, the turbu- 
lent magnetic diffusivity becomes r]i 55 a t c 2 /Ok. 

The sound speed is much less than the Keplerian one is due 
to the disk is cold. It means that h S5 c s /Qk ^ r, i.e. disk is 
geometrically thin. In this model, the disk is considered as an 
infinitely thin, conductive plane. 

The turbulent electro-conductivity <x t = c 2 I{Akt]^ - 
c 2 l(Ana t cJi) is associated with the magnetic diffusivity rj t . The 
surface disk conductivity is A — J cr t dz ~ 2hcr t = c 2 /(2na t c s ), 
and the surface magnetic diffusivity is 



2nA 



a t c s 



In thin accretion disks c s /Vk = h/r <k 1. Thus, an ac- 
ceptable coefficient of the magnetic surface diffusivity is £ = 
(0.01 4- 0.6)(c s /V k )Vk. 

2.3. Initial conditions 

We suggest that, at the initial time, the stellar magnetic field 
with dipole-like topology and magnetic moment p„ penetrates 
the corona and the disk. So, the components of the magnetic 
field B are 



Br 



2p* cos 
~R~i ' 



Bn 



p.* sin 6 



B w = 0. 



We suppose that at the initial time t — 0, the matter of the 
corona and the disk is in mechanical equilibrium with the force- 
free dipole magnetic field, i.e. gravitational force is balanced 
with "centrifugal" force (liquid particle acceleration) and pres- 
sure gradient. 

The momentum equation for the system Q taking into ac- 
count that particles follow circular orbits is 



7 1 

u re, + -Vp = -VO g 



(2) 



where e r is a unit vector with the direction of the cylindrical 
radius r = R sin 8, that on disk surface it is 9 = 7r/2 and r - R. 

This is due to the magnetic pressure B 2 /8n of the dipole- 
like field B decreases like l/R 6 . Also, the most of the simu- 
lation region is occupied by a relatively dense plasma, where 
gas pressure dominates. This dense plasma prevents opening 
of the magnetic field lines. Let us consider that density p is 
function not only of pressure p but also of cylindrical radius 
r. p — p(r,p). Let us denote V{p,r) = 1/p. The momentum 
equation along the z axis (the projection of equation to the 
Z axis) is 

dp d® g 
V<J>,r)-f + — - =0. 
oz oz 

Integrating it by z and suggested that p — > under z — > °° and 
O g — > 0) we finally get that 



p 

I 



V(p',r)dp' + <& g = 0. 



(3) 



The integral is assumed to converge on its interior limit. 

Let us consider now the momentum equation along the 
radial direction 



, dp 30„ 

■ (J 2 r + V(j>,r) 1 f- + — 1 =0, 
or or 



(4) 



After differentiating equation Q by r and subtracting (0} from 
it we get 



r 





orr+ \ ~^ d p' 



0. 



(5) 



Thus, given the function V(p, r), the relation (0 defines the 
function p (r, z) and relation - function oj(r, z). 
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If V(p, r) = k(r)/p a , a = const < 1 (this condition is essen- 
tial for the convergence of the integral in left part of equation 
as p — > 0), then equations Q, l|5} take the form 



f k P 



l-a 

l-a 

k' + u) 2 r = 0. 

{ 1 — a 

Integrating it, we found the function k(r) 



Ink 



GM* J 



n 2 (r)r 2 dr. 



(6) 



Here £2(r) = co(r, 0) is the angular velocity at equatorial plane 
z = 0. 

The angular velocity at the equatorial plane was chosen in 
the following way (rj , r 2 are the parameters of the problem and 

r\ = 2R in , r 2 = 3R in ): 



Q(r) = 



<2 



y jGMJrj 



n - r\ 



<r <r\, 



-(r- n), n < r < r 2 , 



(7) 



, y/GMJr 3 , r > r 2 . 



As a result, the distributions of the pressure p(r, z), density 
p(r, z) and angular velocity uir, z) are obtained 



P(r, z) 



u> = 



l-a GM t 



' l ~ a k'(r) 



!/(!-«) 



. P(r, z) 



P a (r,z) 
k(r) ' 



a - 1 r 



2.4. Boundary conditions 

The surface currents in the disk (at z = 0) lead to a dis- 
continuity in the magnetic field disk-tangential components. 
Mathematically it means that the following condition is ful- 
filled (Landau & Lifshitz, 1982) 



4-71 

nx(B + - B ) = —i, 

c 



(8) 



where B + , B~ are the magnetic field under and over the disk 
respectively, n is the unit normal to the disk and directed down- 
ward in cylindrical coordinates. Since we suggest that MHD- 
flow is symmetric in relation to equatorial plane and consider 
the problem in the upper half-space, then B + = -B = B and 
<03i gives the following 



2tt 

nx B + — i = 0. 

c 



(9) 



Substituting expression for the surface current in the disk i 
which is due to the Ohm's law at the comoving frame i = 
-A(u -V)xB/c into the equation for current, finally we get 



(m - V - £,n) X B = 0. 



(10) 



Assuming that the angular velocity u> and r, (^-components 
of electric field in comoving frame change negligibly in z- 
direction on scales of the order of the disk thickness 2h, we 



obtain for the disk surface electric conductivity A — J cr t dz ~ 
2hcr t . 

In spherical coordinates the tangential components of J10I 

are 



(v- V k )B e -(w-OB„ = 0, 
uB g - (w - £)Br = 0. 



(11) 



The relations (II It . (expressing Ohm's law for surface current in 
the disk), give two boundary conditions in the equatorial plane 
under 6 = n/2 (z = 0). 

The real disk has some vertical (in z-direction) structure 
that is not considered here. It is essential that thermodynamic 
parameters of the disk plasma change in the vertical direction 
turning smoothly to the corona ones. So, the boundary condi- 
tions can be set arbitrarily, based on some physically reason- 
able assumptions. 

We assume that matter flows from the disk to the corona 
at small velocity. It could be expected that plasma leaves the 
disk with the velocity less than the slow magnetosonic one. 
Let us accept that the z-component of the velocity is a frac- 
tion a c of the cusp one in this direction. The cusp velocity is 
not greater than slow magnetosonic one and has the same di- 
rection. The condition fl!2i guarantees that the plasma outflow 
does not exceed the slow magnetosonic one, giving the third 
boundary condition for 6 = n/2 (a = Bj ^\np is the Alfven 
velocity, parameter a c (R) < 1): 



w + a c (R)- 



= 0. 



(12) 



The condition (II 21 guarantees also that from any point of the 
disk five characteristics come out. So, one should set another 
two boundary conditions. 

We assume that matter outflow from the disk does not 
change its interior structure, i.e. the disk has sufficiently large 
mass, energy and angular momentum. We assume that the spe- 
cific entropy of the matter flowing out from the disk does not 
vary with time 



e=7r/2 



= S d (R)- 



To get the equation describing the evolution of the magnetic 
flux function on the disk (at 8 = 7r/2) let us use the induction 
equation 



dBg 
dt 



dREa- 



H 



R dR 



= 0. 



(13) 



Taking into account that on the disk Be = d^/iRdR) and 
Ohm's law is fulfilled, the azimuth component of the surface 
current, + c BrIQuct,) = 0. Integrating Jl 3i . we found 

8^ 

— + (RB R = 0. 
ot 

Let us formulate the boundary conditions on the inner 
bound of the simulation region under R = As for boundary 
conditions on equatorial plane, they are set for reasons of phys- 
ical rationality. The main factor to take into account is on the 
one hand to choose arbitrary the position of the inner bound, 
i.e. the value of the inner radius of the simulation region R m , 
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and, on the other hand, quickly (like R~ 6 ) increasing magnetic 
pressure of stellar dipole-like field. R^ n is chosen to be such that 
in some neighborhood of the inner bound, the star magnetic 
field has dominant influence on the plasma dynamics. In other 
words, in this region Alfven velocity is much more than both 
the gas sound speed and the Keplerian one. On the other hand, 
the choice of a very small R[ n is not reasonable from the com- 
putational point of view, because, into the simulation region, 
in this case includes parts of the magnetosphere where Alfven 
velocity is too large. It leads to an essential decreasing of time 
integration step. Since in the considered model it is assumed 
that magnetic field lines are frozen into the star surface rotating 
with angular velocity Q„, on the inner bound of the simulation 
region (and under it, right on the star surface) plasma moves 
along the rotating magnetic field lines. In the frame of the ro- 
tating star (on inner bound) plasma velocity vector is parallel 
to the magnetic field one. Since transforming to the rotating 
frame B does not change, but u transforms into u - V*e^, where 
V* = £l*,R m sin 9. This boundary condition can be written in the 
following way (R = R{ n ): 



O.lp 0.12 0.15 0.17 20 0.22 0.0(5 0.10 0.14 0.19 0.23 27 

6 



(u - V»e J x B = 0. 



(14) 



Condition (I14t implies also that in the rotating frame, the elec- 
tric field in the inner bound is zero. The conditions (I14> give 
two boundary conditions under R - R[ n . 

On the outer bound, at R — R ouU "free" boundary condi- 
tions are set. Such conditions should not influence on the so- 
lutions inside simulation region preventing poloidal magnetic 
field lines opening. 

On the rotation axis, although it is not a bound, symmetry 
of the flow conditions for this axis are set: 

v = w = 0, B^-Bq- 0, *F = under = 0. 

3. Numerical method and results 

For the numerical integration of the ideal MHD equations 
0, we use Godunov-type conservative high resolution scheme 
(Kulikovskii et al. 119991 Yelenina & Ustyugova 2004 1. To guar- 
antee divergence-free magnetic field we apply the same proce- 
dure as Toth ( 120001 . 

The system of equations ([0 is integrated numerically in the 
region R in < R < R out , < 8 < n/2. We take a non-uniform grid 
in the radial direction, and uniform in the polar angle: Ng = 60, 
Nr = 60. The time step of the integration r is restricted by 
Courant condition. 

To verify the method and results we also used two grids: 
120 x 60 and 240 x 120 and first order numerical method to 
estimate the numerical diffusivity of the scheme. To satisfy the 
goals of this paper, i.e. to investigate the regime of interaction 
between the magnetized star and the disk, it is turned out that 
the grid size presented in the paper is quite enough. 

For mathematical simulation of the interaction between 
a magnetized star and its accretion disk we performed sev- 
eral runs for different surface magnetic diffusivity f . Fig. |2j{7] 
present the results for the following values of parameter if: 0, 
0.001,0.005. 




Fig. 2. (a) Background colour shows the initial distribution of 
the specific entropy S (r, z), solid line presents the magnetic pa- 
rameter B - 1 . (b) Background colour shows the initial dis- 
tribution of the angular velocity u>(r, z), thin lines present the 
magnetic flux function *P(r, z). 



Fig. |3 shows the initial system configuration at time mo- 
ment, t — 0. The background color in Fig.|5^ shows the plasma 
specific entropy distribution S (r, z), thick line corresponds to 
the plasma parameter B = ZnpjB 1 — 1. The background 
color on Fig. I2J3 shows the plasma angular velocity distribu- 
tion u>(r, z), streamlines show magnetic field lines - magnetic 
flux function *¥(r, z). 

The evolution of the coronal magnetic field loops in the 
"star - disk" system depends on the surface magnetic diffusiv- 
ity f . One can pick out the characteristic features which are es- 
sential for this process. For all the cases, poloidal field lines are 
pulled out and reconnected periodically (approximately each 
tenth rotation period of the disk inner bound). Then, a mag- 
netic field toroidal component is generated. After reconnec- 
tion, a plasmoid is formed. It is surrounded by closed poloidal 
magnetic field lines along which the poloidal electric current is 
running. Plasmoid is determined by a strong toroidal magnetic 
field and a low gas pressure. The angular velocity inside the 
plasmoid is different from the corona one. t — 50 is chosen for 
presenting the simulation results (the time is measured in disk 
inner bound rotation period). Up to that moment, several re- 
connections of magnetic field lines took place, and next ejected 
plasmoid moves outwards to the outer bound. Previous series of 
reconnections already led to the opening of the field lines close 
to the rotation axis. We should note that the reconnection of the 
magnetic field lines originates due to the numerical magnetic 
diffusivity. Nevertheless, we suppose that the reconnection and 
plasmoid ejection takes place as well for real magnetic diffu- 
sivity. 

Fig. |5^-f show the distributions of some variables for a 
magnetic diffusivity f = 0. The background colour on Fig. [3ji 
shows the distribution of the poloidal current J p = R sin QB^. 
The magnetic field lines are shown also in this plot by firm 
lines. It is clear from these pictures that the initial configura- 
tion has essentially changed and now it is not dipole-like. In 
the disk differential rotation region (r > 2), field lines come 
out of its surface with large slope angle (more than 30°), fulfill- 
ing the conditions for matter outflow from the disk. A toroidal 
magnetic field is generated as a result of the differential rota- 
tion, implying that there is a poloidal current in the corona, 
forming a double current sheet (background dark colour area 
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Fig. 3. Flow distributions for ( = 0.0 at t — 50. (a) Background colour shows the poloidal current J p . The magnetic field lines 
are shown by firm lines. There is a poloidal current in the corona, forming a double current sheet (shown by background dark 
colour region), (b) Background colour shows the specific entropy S . Streamlines present the matter current lines. Plasmoid 
has hot matter surrounding it from the corona and moves up to the outer boundary, (c) Background colour shows the plasma 
angular velocity co, firm line present the level of the plasma parameter f3 = 1 . In the neighborhood of the star poloidal current 
lines approximately coincide with the poloidal field ones. Fig. (d)-(f) show the distribution of the angular moment transport 
in the system, (d) Background colour shows the magnitude of the angular momentum flux, and streamlines show the direction 
of the system angular momentum transport, (e) Streamlines show the angular momentum flux direction carried by matter. The 
background colour shows the magnitude of this vector, (f) Streamlines present the angular momentum flux direction, carried 
by the magnetic field, its magnitude is shown by colour. There are two regions of intense momentum transport. One of them is 
nearby the star pole where transport takes place due to magnetic stresses (f). The second one is above the disk, in the area of its 
differential rotation where the transport is caused by matter flow (e). 



on Fig. In the neighborhood of the star poloidal current 
lines approximately coincide with the poloidal field ones. 

The background colour on Fig.[3j> shows the distribution of 
the specific entropy S . Streamlines present the matter current 
lines. It is clear that the matter flows into the plasmoid area 
from the star and the disk inner part. Plasmoid has hot mat- 
ter surrounding it from the corona and moves up to the outer 
boundary. 

The background colour on Fig. [5J; shows the plasma an- 
gular velocity distribution, firm line - plasma parameter level 
j6 = 1. It is evident from comparing Fig. [2t and Fig.[3J: that 
the angular velocity is practically constant along the field lines, 
especially nearby the star. It means that in this area the genera- 
tion of the toroidal magnetic field does not take place, i.e. there 
are no poloidal electric currents. The magnetosphere, rotating 
with constant angular velocity, gives angular momentum to the 
plasmoid, twisting the matter inside it. 

Fig-Eft _ EJ show the distributions of several quantities that 
describe the angular moment transport in the system. The an- 
gular momentum conservation equation can be found from the 
continuity equation and is 

dpi 

4r + divL = 0, 
at 



here pi = puRsinO is the angular momentum density. The 
poloidal components of the angular momentum flux density are 

L = R sin 6 L)vu p - j , 

where u p , B p are the plasma poloidal velocity and the poloidal 
magnetic field. The first term in the right part describes angu- 
lar momentum transported by matter, the second one does by 
magnetic field. 

The background colour on Fig. 01 shows the magnitude of 
the angular momentum flux, and streamlines show the direction 
of the system angular momentum transport. There are two ar- 
eas of intense momentum transport. One of them is nearby the 
star pole where transport takes place due to magnetic stresses 
(see Fig. [3f)- The second one is above the disk, in the area of 
its differential rotation where the transport is caused by matter 
flow. These processes are presented on Fig. |3^,f in more de- 
tail. On Fig. streamlines show the angular momentum flux 
direction carried by matter. The background colour shows the 
magnitude of this vector. On Fig.[3J streamlines - angular mo- 
mentum flux direction, carried by the magnetic field, by colour 
- its magnitude. It is clear from the Fig. [3J that angular mo- 
mentum from the star is transported by magnetic field mainly. 
Also there is intense momentum transport inside the plasmoid. 
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Fig. 4. Flow distributions for £ = 0.001 at t — 50. (a) Background colour shows the poloidal current J p . The magnetic field 
lines are shown by firm lines, (b) Background colour shows the specific entropy S . Streamlines present the matter current lines, 
(c) Background colour shows the plasma angular velocity 10, firm line presents the plasma parameter level /3 — 1. Fig. (d)-(f) 
show the distribution of the angular moment transport in the system, (d) Background colour shows the magnitude of the angular 
momentum flux, and streamlines show the direction of the system angular momentum transport, (e) Streamlines show the angular 
momentum flux direction carried by matter. The background colour shows the magnitude of this vector, (f) Streamlines show the 
angular momentum flux direction, carried by the magnetic field, its magnitude is presented by colour. 



The distributions for the same variables at t — 50 are shown 
on Fig. EU5] in case of a finite surface magnetic diffusivity. 
Fig. Eh-EF correspond to £ = 0.001, Fig.Et-§- £ = 0.005. 

The Fig.|6t-c show the influence of the magnetic diffusiv- 
ity on the magnetic field topology. Level lines of the magnetic 
flux function *T(r, z) for all cases (a: ( — 0, b: ( — 0.001, c: 
( = 0.005) are chosen at t — 11. At this time, the first recon- 
nection took place for the case £ = 0. It is clear (see Fig|(J-c) 
that the less the disk conductivity (i.e. the more magnetic diffu- 
sivity), the more the distance from the axis and the disk where 
the plasmoid ejection takes place. In case of imperfect disk 
conductivity, field lines are no longer frozen into the disk and, 
twisted by the star, they are shifted outwards. We can say that 
the more the surface magnetic viscosity the slower the evo- 
lution takes place, and the more seldom plasmoids are formed 
in corona. 

The matter flowing out from the disk and the magnetic field 
both take away the angular momentum from it. The whole an- 
gular momentum flux taken away from the disk per unit time 
is 

L = L m +L f = - J p(rx u)udS + ^- J(r x B)BdS. 

Here, the integration is over the disk surface, dS is the element 
of disk surface directed outwards (from the simulation region). 
The first term L m is the angular momentum flux carried by mat- 
ter, the second one Lf is the angular momentum flux carried by 
magnetic field. 



The time dependence of the angular momentum fluxes car- 
ried by matter L m , by magnetic field Lf and their sum L = 
L m + Lf are shown on Fig.^fr-c for different values of 

For £ = (see Fig.^), the process of angular momentum 
transport from the disk to the corona is quasi-periodic with pe- 
riod equal approximately to ten rotation periods of the disk in- 
ner part. After t„ = 45 there is relaxation of the system accom- 
panied by small oscillations of the angular momentum fluxes. 
During this time there no new plasmoids are generated. After 
t = 65, the reconnection process of the magnetic field lines 
is resumed. It can be seen from the angular momentum fluxes 
oscillations (see Fig.^Jt). 

The finite disk conductivity changes the course of events 
(interaction between magnetic star and accretion disk). As it is 
seen from Fig.^-c, the magnetic lines reconnection lasts until 
some time t* (depending on Q. The more the magnetic diffu- 
sivity f , the later the plasmoids ejection begins and the earlier 
ends. We should note that in case of an imperfect disk conduc- 
tivity the coronal magnetic field evolution is qualitative simi- 
lar to the perfect one. The reconnection of magnetic field lines 
takes place at the time corresponding to the maximum angular 
momentum outflow carried by magnetic field from the disk. In 
contrast to a perfect disk conductivity (£ = 0), the activity of 
the reconnection slows down and the system relaxes to such 
a state that no new plasmoids are generated. The disk angular 
momentum carried by magnetic field gets balanced by angular 
momentum carried by matter. Non-frozen in the disk magnetic 
field lines move along the disk and it lead to the increasing of 
the magnetic flux in the disk. 
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Fig. 5. Flow distributions for £ = 0.005 at t — 50. (a) Background colour shows the poloidal current J p . The magnetic field 
lines are shown by firm lines, (b) Background colour shows the specific entropy S . Streamlines present the matter current lines, 
(c) Background colour shows the plasma angular velocity a), firm line presents the plasma parameter level [3 — 1. Fig. (d)-(f) 
show the distribution of the angular moment transport in the system, (d) Background colour shows the magnitude of the angular 
momentum flux, and streamlines show the direction of the system angular momentum transport, (e) Streamlines show the angular 
momentum flux direction carried by matter. The background colour shows the magnitude of this vector, (f) Streamlines show the 
angular momentum flux direction, carried by the magnetic field, its magnitude is presented by colour. 




Fig. 6. (a)-(c) Influence of the surface magnetic diffusivity f on the magnetic field topology at t = 11 when the first reconnection 
took place for the case f = 0. Level lines of the magnetic flux function *P(r, z) are shown by thin lines (magnetic field lines: on 
(a) for £ = 0.0, on (b) for ( = 0.001, on (c) for £ = 0.005, = 151.6 is shown by solid line. The more the magnetic diffusivity, 
the more the distance from the axis and the disk where the plasmoid ejection takes place. In case of imperfect disk conductivity, 
field lines are no longer frozen into the disk and are shifted outwards. 

the case of imperfectly disk conductivity, magnetic field lines 
are no longer frozen into the disk and, twisted by the star, they 
are shifted outwards. The more surface magnetic diffusivity g, 
the slower the evolution takes place and more seldom plas- 
moids are formed in corona. 

In the case of an imperfectly disk conductivity, the coro- 
nal magnetic field evolution is qualitative similar to the perfect 
one. The reconnection of magnetic field lines takes place at the 
time corresponding to the maximum angular momentum out- 
flow carried by magnetic field from the disk. In contrast to a 
perfect disk conductivity, the activity of the reconnection slows 
down and the system relaxes to such a state that no new plas- 
moids are generated. The disk angular momentum carried by 



4. Summary and conclusions 

The influence of the disk surface magnetic diffusivity on the 
"star-corona-disk" system evolution was studied. For all the 
cases, poloidal field lines are pulled out and are reconnected pe- 
riodically approximately each tenth rotation period of the disk 
inner bound. The magnetic field configuration has essentially 
changed from initially dipole-like one. The magnetic field lines 
come out of its surface with large slope angle, fulfilling the con- 
ditions for matter outflow from the disk, in the disk differential 
rotation region. 

The more the magnetic diffusivity, the more the distance 
from axis and disk where the plasmoid ejection takes place. In 



T.G. Yelenina, G.V. Ustyugova and A.V. Koldoba: Numerical simulations of coronal magnetic field loop evolution 



9 




Fig. 7. Time dependence of the angular momentum fluxes carried by matter L m , by magnetic field Lf and total L = L m + Lf for 
different £ : on (a) for £ = 0.0, on (b) for £ = 0.001, on (c) for ( = 0.005; L m is shown by dashed line, Lf is shown by thin line, 
L is shown by solid line. For £ = (a), the process of angular momentum transport from the disk to the corona is quasi-periodic 
with period equal approximately to ten rotation periods of the disk inner part. Then, after the relaxation the process is resumed. 
For f = 0.001 (b) and f = 0.005 (c), the magnetic lines reconnection lasts until some time. The more the magnetic diffusivity 
the earlier the plasmoids ejection ends. The reconnection of magnetic field lines takes place at the time corresponding to the 
maximum angular momentum outflow carried by magnetic field from the disk. In contrast to a case for £ = 0, the activity of the 
reconnection slows down and the system relaxes to such a state that no new plasmoids are generated. The disk angular momentum 
carried by magnetic field gets balanced by angular momentum carried by matter, in particular, \L\ <*c |L/|. 



magnetic field gets balanced by angular momentum carried by 
matter. Non-frozen in the disk magnetic field lines move along 
the disk and it leads to the increasing of the magnetic flux in 
the disk. 

Summing up we can give following conclusions. 

1. In the case of perfectly conducting disk the evolution of 
the coronal magnetic field leads to quasi-periodic outflow 
of the angular momentum from the disk. The interaction 
process occurs with the reconnection of the magnetic field 
lines and plasmoid ejections. 

2. In the case of imperfectly conducting disk the configuration 
of the magnetic field lines is formed such that the angular 
momentum flux carried by the magnetic field from the disk 
becomes balanced by the flux transported by matter. 
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